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, We propose an experimental model using the Laval nozzle of a sonic analogue of Hawking 

■ radiation. We derive the power spectrum of the outgoing wave emitted from the vicinity of 
Ph ' the sonic horizon instead of the created particle number. Our treatment is based on classical 
^ \ theory, and it should make experiments easier to implement. This experimental feasibility 

' ' . is a great advantage of our model. 

I> ■ 

<N : 

Black hole evaporation, so-called Hawking radiation, is one of the most surprising 
' predictions in theoretical physics. Despite the importance of Hawking radiation, 

. its actual observation in the universe is considered to be very difficult. For this 

\ reason, it seems quite natural to seek experiments in the laboratory that can simulate 

■ this interesting phenomenon. In addition, Unruh showed that there exists a sonic 
analogue of black holes. We speculate that something like Hawking radiation may 

, , be observed in the evaporation of these analogous structures. 

^ ■ Of course, black hole evaporation is a quantum phenomenon. Thus, in the case 

O . of the corresponding laboratory experiments, we should maintain the quantum co- 

herence of the entire system to observe processes of particle creation and changes of 
^ijrj' vacuum caused by the existence of the horizon. With this in mind, several exper- 

• • . iments employing superfluid ^He have been proposed. However, the supersonic 

^ \ motion of this fluid along the walls of the apparatus would cause the collapse of su- 

^ ■ perfluidity and quantum coherence, because vortices may be created near the bound- 

aries. Additionally, even in the an experimental system^-* designed to overcome this 
problem, no evidence of Hawking radiation has yet been reported. 

In consideration of the points raised above, we propose an experiment to observe 
a classical analogue of Hawking radiation, ''^ i.e. the power spectrum of an outgoing 
wave emitted from the vicinity of the horizon. Since we are only concerned with 
classical wave propagation in a fluid flow that has a sonic horizon, this experiment has 
nothing to do with quantum particle creation. Thus, we cannot detect the emission 
of particles corresponding to the evaporation of a black hole. However, Hawking 
radiation possesses another striking feature, a thermal spectrum. It can be shown 
that the power spectrum of an outgoing wave obeys the Planck distribution^^ and 
a role corresponding to that of the particle number in the usual scenario of black 
hole evaporation. Because the behavior we describe does not require a quantum 
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coherence of the system, the experiment we propose is much easier to carry out than 
the previous ones using ^He. 

We first give an outhne of the sonic analogue of the black hole in an ordinary 
fluid. We begin with a perfect fluid, which obeys 



P 



(1) 



where p, p and v are the density, the pressure and the velocity of the fluid. Moreover, 
we assume the case of an adiabatic and rotation-free fluid, implying p = Cp"', 7 = 
Cp/Cv, V X V = (i.e. V = V#), where Cp and are constant pressure and 
constant volume specific heat respectively, and <P is the velocity potential. The 
Euler equation (1) can be integrated, yielding the Bernoulli's equation, 

d'P 1 

^ + -v.v + h{p) = . (2) 



Here, h{p) is the enthalpy and Cg = \/ {dp/dp)ad = \/l p/ P is the local velocity of 
sound. Then, we consider p and # to be perturbated forms written p = pQ{l+ip) and 
= ^0 + where po and ^0 ^ire the density and velocity potential of the background 
fiow. After substituting these into (1) and (2), we obtain the following equations for 
the perturbations: 



dtjj 
ot 



0, 

V(lnpo)- V0 







(3) 
(4) 



The background flow obeys the same equations (1), but with po and vq replacing p 
and V. Eliminating ip in the Eqs. (3) and (4), we derive the equation for the velocity 
potential (j) as 



1 



d 



po \dt 



-Ul + V-^'oH^Ui + ^^o-V --V(poV) 



dt 



1 
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0. 



(5) 



This equation can be interpreted as describing a scalar fleld on the background 
geometry with the metric 



9nv 




- (Cs - Vq ■ Vq) 



■JlJ 



(6) 



where cc is a constant. In order to understand the correspondence between this 
geometry and a black hole, let us consider one-dimensional flow with velocity vq = 
(vo(x),0, 0). In the t-x plane, the metric can be reduced to the form 



1 



Cs{xfdT^ + 



dx 



1 _ vojx)'^ 



(7) 



where we have introduced the new time coordinate r defined by dr = dt + -^^^ . 
This implies that the sonic point, at which vo{x) = Cs{x), may play the role of the 



Letters 



3 



horizon. By the same argument as in the case of the original Hawking radiation, it 
can be shown that this fluid model describes Hawking radiation with a Planckian dis- 
tribution, and we can derive the "surface gravity" which determines the temperature 
of the radiation as 

d{cs - Vq) 



dx 



(8) 

sonic horizon 



Here we propose a model using a Laval nozzle that realizes the situation dis- 
cussed above. Consider an axi-symmetric tube with cross section A{x), where x is 
the coordinate along the tube. When A{x) changes only slightly along x, we can 
regard the flow to be uniform in cross section and one dimensional. For stationary 
background flow along the Laval nozzle, from the Eqs. (1) and (2) we can derive a 
relation between the velocity of the flow and the area of the cross section 

(m^-i)^ = ^, (9) 

where M = vq/cs is the Mach number. This suggests that subsonic flow will be 
accelerated when the cross section becomes smaller, and supersonic flow will be 
accelerated when the cross section increases. In this paper, we consider a tube that 
has narrow throat near the center of the tube and increases at the ends. Such a tube 
is called a "Laval nozzle" . In a Laval nozzle, A{x) decreases along the direction of flow 
{dA < 0) in the region upstream from the throat and increases {dA > 0) in the region 
downstream (Fig.l). Therefore, if we can realize M = 1 in the throat region, we can 
realize supersonic flows in the down stream region and a sonic horizon emerges in the 
throat. Suppose that there exist an asymptotic region 
in the upstream of the flow where the fluid is at rest, 
i.e. Vq = and the pressure and density are given 
by Pq — Pu, po — Pu- It is obvious that this region M<1 
can be related to an asymptotically flat spacetime re- 
gion in the corresponding black hole spacetime. Then 
Bernoulli's equation (2) becomes 

1 2 , 1 Po _ 1 Pu .-.„^ 

2^0 + ^--^- UUj Pig 1 Laval nozzle. 

which determines the velocity of the fluid as a function of the density po/pu- 

Figure 2 displays profiles of the pressure and Mach number along the nozzle 
whose shape is depicted in Fig.l. Under the condition that the inlet pressure at 
X = is fixed, the nature of the flow is determined by the value of the outlet pressure 
at X = 3. We note that a sonic horizon at the throat can be formed without any 
fine-tuning, so that the horizon is stable with respect to changes in the parameter 
values of the nozzle. Furthermore, even in the unsteady situation, experimental and 
numerical studies of time-dependent flow within the nozzle show that the flow is 
settled to the above static one which is predicted by the Bernoulli equation (10). 
Next, let us consider the fluctuations of the fluid. Substituting x) = e"*"^* 
^ K^)d.x -j^^Q j^q^ g^j^j using the WKB approximation, we obtain 

[cs"^ -vl)k'^ -2tovok-to'^ = Q, (11) 
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Fig. 2. Profiles of the pressure p and Mach number M of the background flows along the Laval 
nozzle (Fig.l) in the case of 7 = 7/5 . The inlet, throat and outlet of the nozzle are located at 
a; = 0, 1, and 3, respectively. The vertical dashed-line indicates the shock surface. Above the 
critical value of the outlet pressure (case c), the flow along the nozzle is completely subsonic 
(cases a and b). Contrastingly, the sonic horizon (M = 1 at a; = 1) and regular supersonic flow 
are realized for sufflciently low pressure at the exit (case e). In the intermediate case (case d), 
a shock is formed in the upstream region, which does not affect the flow near the sonic horizon 
at a; = 1. 



whose solutions are given by A; = w/(cs — vq) = fcout and k = —uj/{cs + vo) = 
k\^. The outgoing solution /cout represents a wave that propagates upstream against 
the background flow. In the asymptotic region {vq ~ 0), this mode is a usual 
outgoing plane wave and should be regarded as an observable. However, we note 
that this solution corresponds to a wave that marginally escapes from the sonic 
horizon vq = Cs, so that the behavior near the horizon differs drastically from that 
in the asymptotic region. Contrastingly, /cjn is an ingoing wave that propagates 
downstream and goes through the horizon almost senselessly. Hereafter, we are 
mainly concerned with the outgoing solution k^ut- 

To proceed with the investigation, we assume that the background flow of the 
fluid satisfies p^jpu = g{x/L) and Pq/pu = g^{x/L)^ where L is the characteristic 
length scale of the Laval nozzle. Although the function g is determined by the 
actual shape of the Laval nozzle, its precise form is not necessary in the subsequent 
analysis. We specify only the locations of the throat and the asymptotic region as 
follows: The throat is located at xth (where "th" means "throat"), which satisfies 
g{xtii/L) = (2/(7 + l))^/('^~^), and the asymptotic region in the upstream region is 
given by X ^ +00, fo 0. Under these conditions and using Eq. (10), we obtain 



When we introduce a coordinate z near the sonic point Xth as x = Xth + z , we can 
calculate the surface gravity using Eq.(8), 



^ = ^1^^) ' 5th 



L \ 2 J ' "^^^ dy 



(13) 



Thus, we obtain the following expression for the outgoing mode of the perturbations 
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near the sonic horizon, x ~ Xth 

C"* = exp ( i— In z] = z'^ I (14) 
We note that this wave behaves as a usual out-going plane wave 



lOUt 



Ijj 



exp ( i — X ) (15) 



in the asymptotic region, i.e. x — >■ -^oo. 

In Fig. 3, we give a comparison between the usual black hole spacetime and the 
fluid flow in the Laval nozzle. In the scenario of Hawking radiation, ^) a normal 
mode corresponding to the vacuum, i.e. a positive frequency part, is prepared in 
the past null inflnity 1~ . Here, we consider the evolution of this mode toward the 
past horizon H~ . The existence of the black hole region and the future horizon 
do not affect this evolution so that a quantum state on is nothing more than 
the vacuum state in /~. An observer prepares the vacuum state in the future null 
inflnity which is determined by a positive frequency normal mode in I"*", and 
he deflnes a particle in the vacuum. As in the previous case, we investigate the 
evolution of the mode in backward in time toward the past horizon H~. In this 
case, however, the future horizon signiflcantly affects the evolution of the mode. 
Thus, a comparison between the normal modes in /~ and I"*" in the same place H~ 
shows that the future vacuum in differs from the past one in /~ , and the thermal 
emission of particle can be observed in the future null inflnity I~ . 

We note that the propagation of the perturbation in the background flow in the 
Laval nozzle possesses an essential feature of Hawking radiation in black hole space- 
time. The outgoing plane wave (15) in the asymptotic region, which corresponds to 
the normal mode in the future null inflnity I"*", is affected signiflcantly by the sonic 
horizon, and therefore, its behavior near the sonic horizon is drastically changed, 
as described by (14). It is well known that the logarithmic behavior of the mode 
(14) results in the thermal property of the spectrum of the created particles and its 
temperature is related to the value of the surface gravity k. This means that if we 
could realize the ideal system in which the quantum coherence is being maintained 
perfectly, we would be able to simulate Hawking radiation in the fluid in a Laval 




Fig. 3. Comparison between the black hole spacetime and the flow of the Laval nozzle. 
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nozzle. However, it seems that this would be a very difficult experiment to carry out 
with a sonic horizon, because interaction of the fluid with the wall of the apparatus 
would cause the vortex formation and destroy the quantum coherence of the system. 

Because of the difficulty described above, we treat our problem classically and 
look for the classical analogue of Hawking radiation. ^) We prepare an outgoing wave 
at the past (sonic) horizon H~ as 

•F^" = expiinz) (16) 

which mimics the vacuum state function in the past null infinity . Suppose that we 
observe this wave in the asymptotic region in terms of a power spectrum instead 
of the particle number. For this purpose, we evaluate the Fourier components 

/(a;) = r dz O'i.^r* = r(l - — ) e"^ (17) 

Jo CgK 

of the prepared wave iP'™ with respect to the observed wave in the asymptotic region 
0°"*. Finally, we obtain the power-spectrum that would be observed upstream: 

, „ , , ,9 27ra; 1 

\f(io)r cx —. -. (18) 

This is the Planckian distribution, which can be regarded as the classical analogue 
of Hawking radiation. From Eqs.(13) and (18), the spectrum is characterized by a 
wavelength A ~ 

Of course, our proposal does not represent a full quantum simulation of Hawking 
radiation. However, the nature of the wave propagation near the horizon in the 
presently considered system is similar to that in the case of Hawking radiation. In 
particular, we should be able to observe a Planckian distribution, which is one of 
the important features of Hawking radiation. One of the main advantages of the 
present system is that actually carrying out an experiment and observing the power 
spectrum would be much easier than the previous proposal. In the present case, 
the sonic horizon is believed to be stable. In addition to experimental and numerical 
results^) that indicates this stability, our result (11) itself, which indicates the non- 
existence of growing modes in the perturbation equation (5), can be recognized as 
proof of the stability. Taking into account these advantages, we regard our proposal 
as a first step toward simulation of Hawking radiation in the laboratory. 
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